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Electric generating resources that are part of the interconnected power system must be controlled so that 
their aggregate output matches the electric load at any given time. Frequency deviations and 
unscheduled tie-line transfers between control areas are indicators of a mismatch between load and 
generation. Reliable and economic operation of interconnected power systems relies on multiple layers 
of automatic control systems to insure that generator output follows the changes in load. At the same 
time, these systems must keep the operation of individual system components within safe limits and 
maintain adequate security margins for likely contingencies. 


The variable and mostly uncontrollable nature of wind generation introduces new features into the power 
system control problem. Integration of a continually fluctuating, uncontrolled generation resource such as 
wind impacts the control systems in all time frames. Consequently, the operational and scheduling 
systems must adjust generating patterns to accommodate the variability in the wind in order to maintain 
the same level of system reliability. These adjustments are necessary to ensure that sufficient generation 
is available to meet the control area load and interchange schedules on the various control time frames. 
Among the more significant impacts of integrating large-scale wind plants are sub-optimal generation 
dispatch schedules and the provision of additional generating reserves, both of which affect all control 
system timeframes. The Utility Wind Interest Group (UWIG) has been coordinating a research project to 
quantify these impacts of integrating wind generation. UWIG has secured funding for this work from Xcel 
Energy, Bonneville Power Administration (BPA), National Rural Electric Cooperative Association, 
American Public Power Association, EPRI, Western Area Power Authority, and from its own membership. 
The project is comprised of case studies of the impacts of integrating bulk wind on the Xcel Energy — 
North (Xcel) and BPA systems. This paper summarizes the methodologies developed and initial results 
obtained for the case studies. 
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The operations impact study for Xcel is conducted for 3 different time scales consistent with Xcel’s 
scheduling and real-time control operating procedures. 


= 3-day (72-hour) study horizon of hourly resolution for performing unit commitment simulations 
using operating costs as evaluation criteria 


« 1-hour study horizon of 5-minute resolution for performing intra-hour load following simulations 
using operating costs as evaluation criteria 


« 1-hour study horizon of 4second resolution for performing load frequency control (LFC) 
simulations using statistics of Area Control Error (ACE) as evaluation criteria 
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Xcel’s control area includes 75% of the power consumption of Minnesota and smaller parts of Michigan, 
Wisconsin and North and South Dakota. Xcel’s generating resources are predominantly thermal, with the 
total thermal generating capacities more than 7000 MW. Peaking units fueled by natural gas and oil 
comprise 20% of this capacity. The remaining 80% of Xcel’s thermal capacity is obtained from more 
economical units fueled by coal, nuclear energy, and wood. 


Xcel’s load demand is highest in summer, moderate in winter, and the lowest during spring and fall with 
approximate MW levels of 8000, 6000 and 5000, respectively. Because Xcel’s control area load exceeds 
the capacity of its economic generation, Xcel keeps these economic units online virtually all the time 
except when they are on maintenance or forced outage. Xcel purchases energy to make up for its 
shortage of economical generation to meet its load. During the high-load summer season, the hourly 
import is as high as 2600 MW with a net energy import of about 16% of the total load MWh. Xcel also 
imports energy in other seasons, but the magnitude of the imports is less for these seasons. 


Xcel Energy — North is a member of the Mid-Continent Area Power Pool (MAPP), which provides for the 
sharing of contingency reserves among pool members. Xcel’s reserve share is 161 MW of spinning 
reserve and 160 MW of non-spinning reserve. Xcel also belongs to the Midwest Independent 
Transmission System Operator (MISO), in operation since February 2002 to provide transmission- 
scheduling services. Nonetheless, Xcel currently operates in a regulated, vertically integrated 
environment and is responsible for providing reliable supports in delivering energy to loads within its 
service area. 


Xcel currently has an aggregate wind plant within its control area of about 280 MW capacity located at 
Lake Benton, Minnesota. Annual capacity factor of the wind plant is about 30%, with a seasonal high 
value of 40% for spring and a low value of 15% for summer. Historical data of the wind energy production 
shows a modest diurnal pattern with slightly higher production at night. In the unit commitment phase of 
Xcel’s operations planning, wind generation is discounted by a certain amount from the forecasted value 
to hedge against the inaccuracy in forecasting, especially during volatile wind conditions. Currently, other 
than the real-time generation data, the Xcel control center has no communication link to more accurately 
predict the near-term production of the wind plant. Xcel takes the wind generation on an as-given basis. 
Xcel’s automatic generation control (AGC) reacts to the wind generation variation in real time by 
regulating the ACE. 
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In order to characterize the impact the local wind generation has on system operation, two models related 
to wind-plant operation were developed. The first model was developed for synthesizing the wind 
generation time series over a time horizon with a given time resolution. The second model was developed 
to provide the wind generation fluctuation statistics that determine the reserve requirements used in 
electric utility operations planning. 
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A model was developed for synthesizing the wind generation time series that was statistically 
representative of the actual Xc el wind generation of the Lake Benton wind farm. This statistical model 
was developed from historical wind generation data collected by Xcel as well as high-resolution wind 
generation data collected as part of the DOE/NREL Wind Farm Monitoring Project. The model that was 
developed synthesizes time series for supporting the simulations on the three time scales for the Xcel 
operation study. 


The mechanism of synthesizing the wind generation time series is based on a probabilistic model 
quantifying the probability distribution of wind generation for a time period given its past history. Based on 
the modeling approach adopted by NREL, a Markov chain [1] is used to model the random process of 
wind generation where only the most recent observation is relevant in the probability distribution of future 
wind generation. The decision of using the Markov chain in describing the inter-temporal probabilistic 
behavior of wind generation is based primarily on the simplicity of the model while providing a reasonable 
degree of accuracy, as indicated in NREL publications. Under the Markov assumption, probability 


distributions of the wind generation transition from a given state to different states can be expressed in 
the state transition probability matrix (STM) form: 
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where the probability of wind generation of period t+1, being in state j, given that wind generation of 
period t is in state i is denoted by pi, and each state represents a different MW step for discretizing the 


set of values that the wind generation can assume. Each row of the matrix corresponds to a state at time 
period t, and the row elements are the probabilities of different states at t+1 reached from the state at t. 
Hence the row sum of each row is equal to 1. We further assume the same probability matrix for all time 
periods. 


In order to synthesize wind generation time series for system operations simulations, consider that a 
starting generation level is assigned to i’ at time period t. The row of the state transition probability matrix 
for time period t corresponding to the current state gives the conditional probability distribution of 


generation at t+1; that is p j = 0,1,2... A random number generator of uniform distribution is then 


used to generate a point between 0 and 1. This number is then mapped through the use of conditional 
probability distribution to obtain the state j’ for time period t+1. Repeating the same process with 
generation at t+1, one can then obtain the generation level at t+2 and then up to the end of the study 
horizon for one sample time series. 


In order to validate the probabilistic model, distributions and statistical quantities (mean and standard 
deviation) were derived for both measured and synthetic data and compared. Figure 1 shows the 
probability distributions for 8,760-point measured and 60,000 point-synthesized hourly time series (Note: 
Increased number of synthetic samples have to be produced to simulate all possible combinations of 
state transitions found in the measured data). Table 1 provides the summary statistics for the probability- 
distribution functions shown in Figure 1. This comparison indicates that the developed model produces 
synthesized hourly wind generation time series that are representative of the measured data. Similar 
validation statistics were compared for the higher resolution model, with similar results. 
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Figure 1 Probability Distribution Function for Hourly Measured Data and Synthetic Series 


Table 1. Hourly Standard Deviation and Mean Calculated from Measured Data and Synthetic Series 


Synthetic Measured 


Standard Deviation (MW) 64.0 63.7 
Mean (MW) 71.5 71.6 
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The methodology utilized to characterize the wind generation fluctuation was based on an approach 
presented in [2], in which regulating reserve is considered. The approach is based on using the wind 
generation and system load high-resolution historical data to form the distribution of fluctuation. 
Separately for load and wind generation, the 30-minute moving averages are calculated using historical 
data. Then a time series of 1-minute resolution representing the high-frequency fluctuation is extracted by 
taking the difference between the original historical data and the moving average. Reserve requirements 
are then calculated as 3 times the standard deviation of the high-frequency fluctuation distribution. 
Without wind generation, the requirement is based on the load distribution only. With wind generation in 
place, the combined effect of the uncorrelated fluctuations of load and wind generation results in a 
distribution for which the variance is the sum of the variances of the two components. Applying this 
method to the measured high-resolution wind generation data of April 2001, which is high both in wind 
energy production and fluctuation, it is found that the reserve requirement is increased about 4.8% for 
supporting the integration of current Xcel wind plant of 280 MW into the power grid. 
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The framework for evaluating the operational impacts of wind generation integration encompasses 
simulations on three different time scales: (1) unit commitment time scale, (2) intra-hour load following 
(economic dispatch) and (3) load frequency control (LFC) time scale. Traditional utility scheduling and 
operation tools are used in the evaluation. 


Due to the varying nature of wind, a Monte Carlo approach was adopted. The probabilistic STM-based 
wind plant model tool described in the previous section was used to synthesize multiple wind generation 
time series. The impacts of each wind generation time series are evaluated in a deterministic manner. 
The distribution of the impact values is compiled and used to provide a more representative assessment. 
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Utilities use unit commitment software to schedule generation levels and energy transactions needed to 
meet system load on an hourly basis. Two cost impacts of wind generation integration are identified that 
can be evaluated under this time scale. The first impact is the additional cost associated with inaccurate 
wind forecast estimates for scheduling. Xcel performs unit commitment each day for a 3-day horizon to 
obtain optimal generation and transaction schedules to meet the hourly forecasted load and existing 
transaction schedule profiles. Forecasted wind generation is included in this scheduling process. The 
differential in the forecasted wind generation value used in scheduling and the actual hourly wind 
generation realized requires that in real-time and near real-time operations, other generators be re- 
dispatched from the levels determined in the original schedule. In actuality, the differential will be covered 
in real-time operation by either deploying operating reserves (for deficient wind generation) or by backing 
down economic units (for excess wind generation). Additionally, for a deficiency in wind generation, the 
operating reserves that are deployed in real-time operation must be replaced for subsequent hours, i.e., 
the unrealized wind energy must be replaced from other generating options. Immediate replacement of 
this generation can only be achieved from more costly options (fast starting peaking units, real-time 
markets, etc.) until the real-time operators have time to obtain more economic options (startup additional 
economic units, purchase forward contracts, etc.). Consequently, the generating schedules actually 
implemented in real-time operation are more costly than the schedule determined by the unit commitment 
based on perfect forecasting. This cost difference is assessed by running unit commitment simulations 
for perfectly forecasted wind and inaccurate forecasts. The cost difference between the imperfect and 
perfect forecast runs represents the cost due to forecast inaccuracy. 


The second cost impact that can be evaluated in the unit commitment simulation is the reserve 
component of the load following cost. The generation schedule obtained from the unit commitment 
solution is determined such that the hourly average generation level is sufficient to meet the expected 
hourly average load. However, the control area load is varying continuously in real-time. The utility must 
ensure that sufficient generation is available to cover the sub-hourly changes in load. The load following 
component of these sub-hourly variations is the slow variation associated with the general correlation in 
different customer loads that define the daily load cycle. This variation is on the time scale of several 
minutes, corresponding to the cycle time of economic dispatch execution in utility operation. For example, 
the typical weekday load cycle for Xcel is characterized by a steep ramp-up in the morning from 
approximately 4-10 a.m. and a steep ramp-down in the evening from 8 p.m.-12 midnight, with relatively 
flat load levels throughout the midday and midnight hours. During an hour when the load is ramping 
throughout the hour, the hourly generation level scheduled from the unit commitment will exceed the 
actual intra-hour load value for approximately half of the hour, but will be deficient to meet the intra-hour 
load value for the other half of the hour. Consequently, reserves must be available to deploy within the 
hour to ensure that sufficient generation is available to meet the ramping of the load. In general, the 
amount of reserve required depends on how steeply the load ramps during the hour. The variability of the 
hourly wind generation will affect the amount of reserves that must be available to follow such load 
trends. If the system wind generation follows a daily cycle that is similar to the load cycle, the total load 
following reserves required should be reduced. If, however, the wind follows a pattern that is adversely 
related to the load cycle, as is often the case for strong diurnal wind patterns, the wind would increase the 
intra-hour load following reserve requirement. 


Load following reserve represents the extra generation capability to meet the intra-hour load changes. 
This is an addition to the capacity required to meet the hourly average load. Making such reserve 
available usually requires bringing more units online or scheduling unit generation levels in a less 
economic manner in the operation planning stage. This results in extra cost. 
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The deployment of the available load following reserve to meet the intra-hour slow variation of load 
changes also results in extra cost. We refer to this cost as the energy component of the load following 
cost. An economic dispatch program is used to simulate the intra-hour deployment of generation every 5 
minutes for 1 hour. The unit on/off statuses are fixed according to the solution of the unit commitment run. 
The economic dispatch simulation models the contingency reserve requirement and regulating reserve 
requirement. It models the load following reserve requirement dynamically in the sense that at a given 
time step, the amount of reserve is reduced to match the increase in load or reduction in wind generation. 
The economic dispatch is in essence deploying the reserve and converting it into generation. The 
economic dispatch program used for this study also models an artificial unit that is dispatched when load 
and reserve requirements are not met by the actual units currently online. A penalty charge 
commensurate with the peaking unit average generation $/MWh cost at full load is assessed for the 
dispatch of this artificial unit, representing an equivalent use of peaking energy or purchase of spot 
market energy to meet the load requirement. 
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In addition to providing generation to meet the slow intra-hour variations associated with the ramping of 
the load through the load cycle, the minute-to-minute variations in load must also be accommodated. 

This regulation service requires that fast-responding generation be reserved to respond to these fast 
fluctuations. The minute-to-minute variability of the wind production will impact the regulating reserve 
amount required to meet these changes. As discussed in the previous discussion of the Reserve 
Requirement Model, it was determined that the additional amount of regulating reserve required to 
support the currently existing wind generation penetration of Xcel is minimal. Therefore, it is assumed that 
no additional regulating reserve is required. Simulations of the real-time regulation process with wind and 
without wind were performed, however, and the control performances determined by the simulations were 
compared. 
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Figure 2 provides a graphical overview of the processes utilized to evaluate the impacts of wind on NSP 
system operations for the three separate time scales discussed above. The upper-layer of the flow chart 
shows the use of unit commitment to determine the wind generation forecast inaccuracy cost. A three-day 
period was selected for each of the winter and summer seasons where winter is characterized by high- 
wind and medium load level and summer by low wind and high load level. Unit commitment simulations 
were performed for 100 corresponding wind generation time series synthesized from the wind STM tool to 
represent the actual wind generation. Additional sets of simulations were performed to determine the 
cost of actual generation when the wind generation forecast is off by 10%, +20%, and +50% of the 
actual wind generation. Comparison of the cost distributions for the original set of simulations relative to 


the cost distributions for the various sets of inaccuracy simulations provide the basis for determining this 
cost component. 


The middle-layer of the flow-chart shows the process for assessing both the reserve and energy 
components of the cost impact of wind on intra-hour load following. Load following reserve requirements 
(LFRR) with and without wind generation are estimated from historical data on the basis of the hourly load 
and wind generation changes. Unit commitment is then used to determine the cost differential under 
different reserve requirements between with and without wind generation scenarios. 


The simulation process for calculating the energy component of load following cost is shown in the lower 
half of the middle layer. The economic dispatch tool is used to simulate the intra-hour effects of re- 
dispatching generator set points every 5 minutes. Rather than simulate all 24 hours for each of the 3-day 
UC periods, the economic dispatch is performed for 4 representative hours from the Xcel daily load cycle 
— hours 3 (nightly flat), 8 (morning ramp-up), 14 (afternoon flat), and 23 (night ramp-down). The 
commitment schedule for these hours is taken from the median wind generation simulation for each of the 
2 wind “season” scenarios. Xcel provided 5-minute resolution load data for 5 days during summer 2000. 
This data provided hourly load shapes that were scaled for the hourly average load values for the 
selected hours from the 3-day simulation windows. For each of the selected hours for each wind season, 
100 wind series of 5-minute resolution and 1-hour horizon were synthesized using the STM wind tool. 
The economic dispatch was performed for each of these combinations of load and wind series for the 
selected hours for the 2 wind seasons to provide distributions of the energy component cost for the 
selected hours. These cost distributions were then used to determine an annualized energy cost 
component. 


The bottom layer of the flow chart shows the simulation process for assessing the impact of wind on 
regulating minute-to-minute fluctuations in load. To perform this assessment, two high-resolution (4- 
second) wind generation series of 1-hour duration were selected from the NREL wind generation data set 
for each of the selected hours for the relevant wind season. Xcel provided 4-second resolution load data 
for the same 5 summer days for which the 5-minute resolution data was provided. For each of the 
selected hours, four high-resolution load series were selected from the Xcel historical data. The LFC 
simulations determined the ACE statistics for each combination of load and wind, which provided 8 
distributions of 900 calculated ACE values. These ACE statistics were also averaged up to 1-minute ACE 
values to correspond to the NERC performance metric period. These statistics were compared to the 
ACE statistics calculated for the load series without wind generation. The results showed that the 
addition of wind had very little negative impact on system ACE. Consequently, it was assumed that the 
cost impact was minimal as well. 
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Figure 2. Flow Chart of the Simulation Framework Utilized for the NSP Wind Impacts Case Study 
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All results presented in this paper are based on the developed models of NSP’s operation planning and 
controls systems and operating procedures for the current system. These models include many 
approximations and extrapolations such as annualizing based on small discrete time periods, which are 
not discussed at length due to paper length requirements. The authors believe that most of these 
approximations result in conservative impact (more significant impact) calculations. The results 
presented are, however, based solely on NSP’s current system configuration and wind penetration level 
of 280 MW on an 8000 MW peak system. At some higher penetration of wind, the load following and 
regulation reserve requirement impacts that were found to be negligible may result in additional operating 
cost. Those sensitivities have not yet been determined, but should be noted when considering the results 
presented. 
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For both winter and summer scenarios, simulations are performed to determine the increase in total 
operating cost for forecast errors of 10%, +20% and +50% of the realized wind generation, with a 
positive error representing an over-optimistic forecast and a negative error an over-pessimistic forecast. 


In evaluating this cost increase, it is assumed that the forecast amount is used to reduce the load in day- 
ahead operation planning for scheduling other generating resources. The following assumptions were 
made for modeling Xcel’s hour-ahead and real-time re-scheduling procedures: 


e Optimistic forecasting of wind generation results in purchasing less energy than necessary in the 
day-ahead planning. Interchange is fixed at the forward schedule with no further transaction being 
arranged. Peaking units are dispatched to compensate for the unrealized wind generation at real- 
time. 


e Pessimistic forecasting of wind generation results in purchasing more energy than necessary in 
the day-ahead planning. Generations from economic units are backed off for the unexpected 
increase in wind generation at real-time. 


Figure 3 and Figure 4 plot the Monte Carlo expected value of forecast inaccuracy cost in % of the total 
cost as a function of the % forecast inaccuracy for the winter and summer scenarios, respectively. The 
winter scenario results for 10% and 20% pessimistic forecasts are small negative values where positive 
values are expected. Such discrepancy likely arises from the finite numerical precision of the unit 
commitment program. 


The over-optimistic forecast inaccuracy results in an increase of generation levels during re-scheduling. 
The increase in total cost for such forecast error relative to the perfect forecast case implies that marginal 
cost for generation increase is strictly greater than the highest transaction price being exercised in the 
day-ahead operation planning. The over-pessimistic forecast inaccuracy results in a reduction of 
generation levels during re-scheduling. The associated cost increase relative to the perfect forecast case 
implies that the marginal cost for generation decrease is strictly less than the highest transaction price 
being exercised in the day-ahead operation planning. Note that the marginal cost for generation increase 
is greater than or equal to the marginal cost of generation decrease. 
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Figure 3. Plot of Expected Total Cost Increase versus Forecast Inaccuracy, Winter Case 
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Figure 4. Plot of Expected Cost versus Forecast Inaccuracy, Summer Case 


The cost increase due to forecast inaccuracy for the summer case is much more pronounced than for 
winter. In the summer case, the marginal cost of the transaction is much higher because of the higher 
load demand. All the inexpensive generating units are scheduled up to their high limit. As a result, for an 


over-optimistic forecast error as small as 10%, a peaking unit must be dispatched for additional 
generation in re-scheduling as demonstrated in the simulation. The sharp increase in cost from perfect 
forecasting to 10% over-optimistic in forecasting is due to the high no load cost and startup cost of the 
peaking unit. 


In real life, forecasting error is random, fluctuating between positive and negative values of varying 
magnitude. The error associated with a given forecast is unknown until after the actual wind generation is 
observed. In evaluating the effect of randomness in forecasting, the forecast error was modeled as a 
random variable with a given probability distribution. Furthermore, the extra operating cost due to error in 
forecasting was treated as a random variable for which its expected value can be evaluated. For 
simplicity in calculation, a highly idealized uniform distribution with equal range for positive and negative 
error values was considered. 


The numerical results of the extra operating costs for different percentages of over-optimistic and over- 
pessimistic presented previously are referred to as inaccuracy cost functions. Under randomness in 
forecast error, the expected extra cost can be calculated by integrating the two inaccuracy cost functions 
appropriately weighted by the probability distribution. Consider using a scaled forecast for unit 
commitment operation planning. The value of expected extra cost is calculated differently as the 
inaccuracy cost functions are weighted by the distribution differently in integration. In fact, it can be shown 
that the expected extra cost can be reduced with respect to no scaling by using some scaling factor less 
than 1, unless the inaccuracy cost for over-pessimistic is substantially higher. For different distribution 
ranges, the following table gives the highest scaling factors such that only the over-pessimistic cost 
function is used in the integration evaluation. Such scaling factors are quite close to the optimal factors, 
which yield the lowest expected extra costs. 


Table 2. Sub-optimal Scaling Factors 


The expected extra costs are calculated for winter and summer using the sub-optimal scaling factors just 
proposed in Table 2. Assuming the other two seasons are about the same as winter, the annualized extra 
costs in $/MWh over all four seasons are shown in Table 3. The extra costs in ($/MWh) are based on the 
MWh values of wind generation time series for Monte Carlo simulations. 


Table 3. Expected Extra Costs for Different Distribution Ranges, Annualized 


Extra Cost (§/MWh) 


Using a scaling factor less than 1 is consistent with the operation planning strategy of Xcel. A scaling 
factor less than 1 means that the wind energy value used in operation planning is reduced from the 
forecast. This is equivalent to using the actual forecast in operation planning, and then setting aside extra 
reserve to hedge against the forecasted wind energy not being realized in real time. 
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Calculation of the LFRR from the historical data indicated that the addition of wind does not significantly 
increase the reserve requirement. Consequently, the reserve component of the load following cost is 
assumed to be zero. Again, it should be noted that this determination is for the existing NSP wind 
penetration level. Assuming a reserve component cost of zero for wind means that the energy 
component assessed using intra-hour economic dispatch simulation will be higher than the energy 
component cost that would be calculated if additional load following reserves were added to support the 
wind. 


As described in the Analysis Framework section, four different hours of the day were selected for 
simulation, representing different load ramping and wind variation characteristics. Results of the energy 


component of load following cost in ¢/MWh, for each of the simulated hours for winter and for summer are 
presented in Table 4. 


Table 4. Load Following Cost for 4 Different Hours of Day for Winter and Summer 


52.289 736.105 


In projecting the average cost over a day, we divide a day into 4 different periods based on the load 
ramping characteristic with each period including a simulated hour. The weighting factors for the 4 hours 
are 2, 8, 8 and 6, respectively. Under such distribution, the load following costs for a typical day for winter 
and summer are calculated by weighted averaging. Assuming that load following cost for spring and fall is 
the same as winter, the annualized cost is calculated by averaging over all four seasons. Table 5 shows 
the cost for winter, summer and annualized. 


Table 5. Load Following Cost for Winter, Summer and Annualized 


Load Following Cost (¢/MWh) 39.38 45.03 40.785 
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Selected hours for the regulation simulation are the same as those for the intra-hour load following 
simulation. 


Table 6 and Table 7 show the simulation results of standard deviations of ACE, which reflect the random 
fluctuation of the variables being controlled. Among the four simulated hours, standard deviations are of 
the same order, though somewhat larger for hour 8. The standard deviations are almost unchanged 
between the without and with wind generation scenarios. This suggests that NSP’s current wind 
penetration of 280 MW on an 8000 MW peak system has no impact on the control performance. This 
means that for NSP’s current wind penetration level and regulating capacity and for the reserves 
allocated in the simulations, the variability of the wind on the 4second time frame didn't significantly 
exceed the capability of the system to follow these variations. Accordingly, the cost impact of additional 
regulating reserves to accommodate wind is assumed negligible. 


Table 6. Standard Deviation of 4-second ACE from without and with Wind Simulation 


Standard Deviation of 4-second ACE 


802 | 8.89072 


Standard Deviation of 1-minute-average ACE 
Hour of Day : - : ; 
Without Wind Gen With Wind Gen 
12.3548 12.3369 
15.3375 15.3297 
10.9075 10.9025 
23 13.3211 13.3477 


Using the results obtained for the 50% forecast inaccuracy range, the total cost impact for the current 
NSP system is found to be approximately $1.85/MWh. 


TER OD SOO UNHAV EVE ea KOH 
AVEO Hee HRD 


Over 1200 MW of wind power is expected to be installed in the Bonneville Power Administration service 
territory over the next few years. The total anticipated wind capacity to be developed would be spread 
over several to many individual projects, with the latest information projecting a total in excess of 1200 
MW in 10 different projects. The majority of the plants are to be located in the mid-Columbia River region, 
south of the John Day “cutplane,” reducing concerns related to existing BPA transmission constraints. As 
part of the UWIG research project for which the Xcel case study was performed, an additional case study 
is being performed to determine how such a significant quantity of wind generation can influence the 
operation of the BPA system, especially the hourly and daily scheduling and coordination of hydropower 
facilities. 


Daily scheduling of BPA hydro resources involves navigating a myriad of non-electrical constraints. 
Violations of any of these constraints have major consequences. The BPA transmission system and 
hydro resources are required to meet the varying demand for electrical power efficiently and 
economically, while continually complying with non-power requirements. The uncertainty associated with 
hourly production of substantial wind resources over the course of a day as well as continuous 
fluctuations within the hour, and the possible impacts on hydro resource operation is the major question 
to be addressed in this case study. 


The operations impact study will encompass three scheduling time frames: monthly, pre-schedule, and 
real-time (next hour). Representing the proposed wind generation on each of these time scales, both on 
an individual plant and an aggregate basis, is a critical element of the case study. Because historical 
data from large windplants in the Northwest is limited, a variety of data sources and techniques must be 
employed to build windplant models for the various calculations and simulations to be performed in later 
tasks. Because of the high-resolution wind speed data limitations and the fact that most of the wind 
plants have not yet been developed, the initial impact study will utilize hourly resolution simulations to 
determine an estimate of the impact once full build out is achieved. This paper describes the developed 
hourly resolution wind plant models. Although the initial impact study is not yet complete, this paper also 
presents the methodology to be implemented for the initial system impacts evaluation. 


HE kc LEE 


Unlike the Xcel wind plant, a very limited data record for characterizing the BPA wind plant currently 
exists. Only 4 of the 10 proposed wind projects that are to be considered in the impact study are online, 
and three of those only since 1/1/2002. The other 6 projects are still in the development stage. 
Therefore, historical wind generation data at these sites is presently unavailable. A method for producing 
expected generation levels for each wind project based upon local wind profiles was developed using 
wind speed data from various points within the region along with wind project-specific information. 


The primary source of wind speed data was the Oregon State University, Pacific Northwest Wind Data [8]. 
This database contains a historical record of hourly-resolution wind speed measurements at many sites 
throughout the region. None of the OSU measurement sites exactly correspond to the BPA wind project 
sites. Based on type of data available, range of historical data, and proximity to the wind project sites, 
three OSU data collection sites were found to provide an adequate representation of the wind profiles at 
the 10 BPA wind project sites. These 3 sites provided synchronized hourly wind speed and direction data 
from 10/78 — 12/00. 


Each of the wind project locations was mapped to one of the three OSU wind measurement locations. 
Based upon the mapping that was utilized, the average distance between a wind measurement site and 
respective wind project site is approximately 18 miles. In order to represent the geographic diversity of the 


proposed BPA wind projects, a technique for transposing the OSU wind speed measurement time series 
from one site to another was developed. The basic concept is that the individual project outputs are 
calculated based on a time-shifted wind speed series. The time shift in each project seeing the wind time 
series is calculated from the relative proximity of the wind project site and the OSU anemometer location 
that recorded the wind speeds. Based on these spatial vectors, the time shift is calculated based on a 
constant wind direction and wind speed at which the series will move from the anemometer location to the 
project site. The wind mass, represented by the wind speed time series, is assumed to move as a 
continuous line orthogonal to the direction of the wind. The constant wind speed and direction are 
determined for each time series as the average wind speed and direction over that period of time. For 
time series longer than 12 hours, the time series is broken into sub-series of 12-hour spans, and the 
transfer procedure performed on each sub-series. 


Using the transposed wind speed measurements, the wind speeds are scaled according to project- 
specific hub height and converted to individual plant production data using plant-specific turbine power 
curves. Based upon investigation of various wind turbine modeling issues, for resolutions greater than 
several seconds, a quasi steady-state model incorporating the wind turbine power curve is sufficient for 
developing wind plant production data. 


For each of the wind farm projects, a scaling factor has been introduced that accounts for various 
conditions that can result in higher synthesized output production than would normally be expected and 
that can not be explicitly modeled, including 


= Wind farm losses that occur within the plant (generator step-up transformer, collection feeders, 
substation transformer) 


« Wind turbines that are offline due to maintenance or repair. Because turbine availability 
information is not available, this cannot be included explicitly. 


Comparisons of measured and calculated hourly plant production data were made with the four projects 
that have been brought online. For these four sites, statistical analysis was performed on actual 
measured BPA wind production data and the calculated wind production using the OSU wind speed data 
sets for the period of 2/1/02 — 3/31/02. As expected, without using the scaling factors mentioned above, 
the mean calculated project outputs are found to be higher than those measured. Therefore in order to 
better match the calculated data with the measured data, a scaling factor was used for each of the 4 wind 
project sites. The scaling factors were chosen such that the mean of the distribution of calculated 
generation of each windplant is approximately the same as the mean of the distribution of measured 
generation values. For the data sets analyzed, the stated scaling factors yield distributions that are more 
similar in shape to the measured distributions and yield wind generation time series that are more similar 
to the actual wind generation profile for the time period. The trade-off with using the scaling factor to 
better match the mean production level is that the calculated wind farm does not peak as high as 
measured data, meaning that the full range of wind generation fluctuations is not represented by the 
model. 


Sample time series plots for a single project site comparing both the calculated and measured project 
power outputs are shown in Figure 5. The time period shown was selected somewhat arbitrarily. The 
calculated power during certain portions of the time series tracks the measured power output very well. 
Also, without the scaling factors mentioned previously, the calculated power would tend to exceed the 
measured by a significant amount. Some of the hour-to-hour fluctuations are also not found in the 
calculated outputs. 


— Calculated 


Figure 5 Sample Measured and Calculated Plant Production Time Series 


Sample distributions for the same project site for both the measured and calculated series are found in 
Figure 6. A visual comparison of the calculated and measured distributions for each of the sites shows 
that the general shape of the distributions for each project is similar. The notable differences in the 
distributions are at the extremities, some of which is a result of the introduction of the scaling factor. Also 
contributing to this fact is that the calculated generation levels tend to exhibit the cut-out at wind speed 
greater than rated when the measured data often doesn’t. The higher percentage found at the highest 
level of calculated power is the result of the wind speeds reaching rated values, and the entire plant 
project outputting rated value (scaled with the scaling factor, however). 


Calculated 
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Figure 6 Sample Measured and Calculated Wind Plant Production Distributions 
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The initial impacts study is still in progress, and as of yet, there are no results to present. The 
development of the framework for the study is near completion. The focus of the initial impact study will 
be to determine the impact of the proposed BPA wind plant on BPA’s hydro operations and on the ability 
to export energy out of the BPA region. The flow-chart shown in Figure 7 graphically summarizes the 
process to be utilized for conducting the BPA preliminary impacts study. Note that the tools to be utilized 
are not the standard operational planning tools that were utilized for the Xcel case study, but rather BPA 
custom tools designed specifically for the requirements of BPA’s hydro-generation system. 
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Figure 7. Flowchart summarizing the analytical framework to be exercised for the BPA preliminary impacts study. 


The general process is as follows. An aggregate BPA wind plant production time series for specific 
“water seasons” will be obtained from the wind model tools described previously. This wind generation 
will be used to reduce a forecasted load series for the same time period to produce a net load series that 
must be served from the BPA system. This net load series and an estimated series of water available to 
the BPA hydro system serve as inputs to the BPA hydro simulation tools. BPA’s 30-day scheduling tool 
will be utilized to produce a hydro generation schedule to meet the net load series while maintaining all of 
the non-power constraints of the hydro system. The hourly hydro schedule produced for each wind 
scenario will be compared with a base case “no wind” hydro schedule to determine the cost impacts of re- 
dispatching the hydro system due to the wind fluctuation. 


The second cost impact to be evaluated in the preliminary impacts study is the extent to which wind limits 
BPA’s ability to export energy during high load seasons. For each simulated hydro generation schedule, 
generation snapshot for specific hours of the schedules will be modeled in a power flow simulation to 
determine the impacts of the additional wind generation on the ability of BPA to export and import energy. 
Any impacts on the transmission capacity associated with the added wind generation will be valued. The 
aggregate cost impacts of wind on the hydro and transmission systems will be compared with no wind 
scenarios to provide the preliminary cost impact estimates. 


“LO ALA 


This paper reports two case studies on the operational impacts of integrating large wind plants into the 
power grid. Although the objectives of the two case studies are similar, the methodologies utilized to 
achieve those objectives differ significantly according to the differences in the environments in which the 
two utilities operate. Xcel Energy —North is predominantly a thermal generation system with an existing 
wind plant of little geographic diversity for which exists a significant historical record of wind and energy 
production. BPA, however, is predominantly a hydro system to which a very geographically diverse wind 
plant is expected to be added. Additionally, there is very little historical record of the wind energy for the 
existing BPA wind plant, and a historical record of wind speed exists only for regional sites not specific to 
the proposed wind plants themselves. 


The general analysis approach for both case studies is to use the tools utilized by the respective utilities 
to simulate the respective operational procedures to determine the cost impacts of integrating their 
respective wind plants. Accordingly, it should be noted that the results of the impact studies presented in 
this paper are specific to the systems studied as they currently exist. Furthermore, it is very difficult to 
exactly model all of the operational scheduling and real-time operation procedures for a given utility and 
to obtain all of the necessary data, so assumptions and extrapolations are made in the developed 
models. The investigators attempted to select these simplifying measures such that the effect was to 
produce a more conservative (more significant) impact. The sensitivity of the results to many critical 
parameters such as wind penetration level, generation mix, and energy transaction pricing should be 
investigated to more fully understand the impacts of integrating bulk wind generation. 
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